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Abstract 

We generalize the fermion Chern-Simons theory for the Fractional Hall Ef- 
fect (FQHE) which we developed before, to the case of bilayer systems. We 
study the complete dynamic response of these systems and predict the exper- 
imentally accessible optical properties. In general, for the so called (m, m, n) 
states, we find that the spectrum of collective excitations has a gap, and the 
wave function has the Jastrow-Slater form, with the exponents determined by 
the coefficients m, and n. We also find that the (m, m, m) states, i. e. , those 
states whose filling fraction is ^, have a gapless mode which may be related 
with the spontaneous appearance of the interlayer coherence. Our results also 
indicate that the gapless mode makes a contribution to the wave function of 
the (m, m, m) states analogous to the phonon contribution to the wave func- 
tion of superfiuid He4. We calculate the Hall conductance, and the charge 
and statistics of the quasiparticles. We also present an SU{2) generalization 
of this theory relevant to spin unpolarized or partially polarized single layers. 

PACS numbers: 7.10.Bg,71.55.Jv 
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I. INTRODUCTION 



The Quantum Hall Effect is observed in two dimensional electron systems (2DES) in the 
presence of strong perpendicular magnetic fields, at very low temperatures. This effect is 
characterized by the existence of an energy gap between the ground state and the lowest 
excited state. In the case of the Integer Quantum Hall Effect (IQHE) the energy gap is 
the Landau level spacing produced by the external magnetic field at integer filling factors. 
In the Fractional Quantum Hall Effect (FQHE) the energy gap appears as a result of the 
interparticle correlations due to the strong interactions between the electrons. 

If one allows for the presence of new degrees of freedom, a richer variety of states can 
be found. The two obvious possibilities that one can consider are systems in which the 
electronic spin is not frozen by the Zeeman energy, and systems in which two or more layers 
of 2DES are coupled together. For instance, the experimentally observed = | state 
has been explained theoretically by Haldane and Rezayi using the fact that the system 
is not spin-polarized. 

Due to continuing advances in material-growth techniques, it has been possible to fab- 
ricate high-quality multiple 2DE layers in close proximity. In these systems the layer index 
is the new degree of freedom, and the interplay between the intralayer and the interlayer 
Coulomb interactions gives rise to very interesting physics. In particular, this competition 
can explain P] the experimental observation [Q of the destruction or weakening of the IQHE 
at odd filling fractions. Another interesting case is the one of the = | state. In single-layer 
systems, even though many transport anomalies have been reported, there is no evidence of 
FQHE. On the other hand, this is a well observed FQHE state in double-layer systems. 

Motivated by the fact that very interesting physics can be found in these 2DES if one 
considers new degrees of freedom, we study double-layer FQHE systems. Our formalism can 
also be extended to the study of spin non-polarized systems. 

There are two energy scales that play a very important role in this problem. One is the 
potential energy between the electrons in different layers, and the other one is the tunneling 
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amplitude between layers. We only consider the case in which the tunneling between the 
layers may be neglected, and both layers are identical. Therefore, the number of particles 
in each layer is conserved, and the collective modes corresponding to in phase and out of 
phase density oscillations are decoupled. 

We generalize the fermionic Chern-Simons field theory developed in reference [Q. The 
generalization is straightforward. We consider a theory in which the electrons are coupled to 
both the electromagnetic field, and to the Chern-Simons gauge fields (two in this case, one 
for each layer). We show that this theory is equivalent to the standard system in which the 
Chern-Simons fields are absent, provided that the coefficient of the Chern-Simons action is 
such that the electrons are attached to an even number of fluxes of the gauge field in their 
own layer, and to an arbitrary number of fluxes of the gauge field in the opposite layer. In 
this form, the theory has a U{1) ® U{1) gauge invariance. We obtain the same action as the 
one derived by Wen and Zee in their matrix formulation of topological fluids [0. 

In this paper, we study the liquid-like solution of the semiclassical approximation to this 
theory. We can describe a large class of states which are characterized by filling fractions in 
each layer given by 

n - (±^ + 2.2) 

z/i - 



n2 - (±^ + 2si)(±^ + 2S2) 



' n2-(±^ + 2.0(±^ + 2.2) ^'-'^ 
where pi, p2, si, S2, and n are integers. This includes the so called {mi, 7712, n) with filling 
fractions u = 2n-mi-m2 ^-j^^ (m,m,m) states, with filling fractions u = —. We calcu- 

ji'^— mim2 ' \ 7 7 / 7 o ^ 

late the electromagnetic response functions and find the spectrum of collective excitations. 
We find that for the {mi, 1712, n) states the in phase as well as the out of phase collective 
excitations are gapped. On the other hand, for the (m, m, m) states there is an out of phase 
gapless mode which, in the absence of tunneling between layers, indicates the spontaneous 
breaking of the U(l) symmetry associated with the conservation of the relative number of 
particles. 
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We show that, already at the semi classical level of our approach, the density correlation 
functions saturate the /-sum rules, associated with the two separated conservation laws, 
i. e. , the number of particles on each plane is separately conserved. Using this property, 
we can derive the universal form of the absolute value squared of the ground state wave 
function at long distances and in the thermodynamic limit. For the (m, m, n) states the wave 
function that we find has the Jastrow form predicted by Halperin [§. For the (m,m, m) 
states, we find that there is an additional factor which represents the oscillations of the 
gapless mode. This additional factor has the same form as the contribution of the phonons 
to the superfiuid He4 wave function. Exactly as in the superfiuid He4, the gapless mode 
factor gives a negligible small contribution to the ground state energy but it is crucial to 
get the correct correlations. In fact, MacDonald and Zhang have recently calculated the 
collective excitations for double-layer FQHE systems using the single mode approximation 
on the Halperin (m, m, m) state , and find that it violates the sum rules. 

We calculate the quantum numbers of the quasiparticles. We find that the charge of 
the quasiparticles for a given layer, is determined by the filling fraction of the layer, and 
by the effective number of Landau levels filled in that layer, at the mean field level of the 
semiclassical approximation. We compute the statistics of the quasiparticles. We show that 
it is well define only in the case in which the coefficient of the Chern-Simons term in the 
effective action for the fluctuations of the Chern-Simons gauge fields, is not singular. In 
this case, the statistics of the quasiparticles of a given layer is proportional to the number 
of fiuxes attached to the electrons in that layer. On the other hand, the relative statistics 
between quasiparticles in different layers, is proportional to the relative number of fluxes. In 
the case in which the statistics is not well defined, as for instance for the (m, m, m) states, 
we find that the effect of the gapless mode is to induce a long-range attraction that forces 
the quasiparticles in different layers to move together, forming a bound state. 

We also develop a generalization of our theory in which the SU (2) spin symmetry is 
taken into account explicitly. We obtain the general form of the SU (2) hierarchies and show 
that, in general, they do not coincide with those of the U{1) ® U{1) theory. 
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The paper is organized as follows. In Section |T| we discuss the generalization of the 
fermionic Chern-Simons field theory to double-layer systems, and derive the electromagnetic 
response functions. In Section ^ we calculate the spectrum of collective excitations for the 
(3, 3, 1) state and for the (m, m, m) states. In Section |V| we show that the density response 
functions calculated within the gaussian approximation saturate the /-sum rule. We also 
derive the form of the absolute value squared of the ground state wave function for these 
states, valid at long distances and in the thermodynamic limit. In Section we derive the 
Hall conductance and show that it has the correct value already at the semiclassical level of 
our approximation. In Section |V| we discuss the quantum numbers of the the quasiparticles. 
In Section [VII| we discuss the SU (2) version of this theory. Finally, in Section |VIII| we 
summarize our results. Two appendices are devoted to the proof of fermion-to-fermion 
Chern-Simons mapping in bilayers and to the estimate of the contribution of the gaussian 
fluctuations of the ground state energy of the (1, 1, 1) state. 



II. FERMIONIC CHERN-SIMONS THEORY FOR DOUBLE LAYER FQHE 

SYSTEMS 

In this section we describe the generalization of the Chern-Simons field theory for the 
single-layer FQHE to a double- layer two dimensional electron system (2DES). 

In the second quantized language, the action for a double-layer 2DES in the presence of 
an external uniform magnetic field B perpendicular to it is given by 

S = jd'zY: {C(^)[^^0 - l^W.iz)?] 

-\jd?Z J d?z' Y: {\Mz)\' - Pa)V^fs{\z- z\){\M^')\' - Pf^) (2.1) 

where the indices a = 1,2 and (3 = 1,2 label the layers, pa is the average particle density in 
the layer a, ip{z)a is a second quantized Fermi field, /Iq is the chemical potential and is 
the covariant derivative which couples the fermions to the external electromagnetic field A^. 
In what follows we will assume that the pair potential has either the Coulomb form, i.e., 
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vmu = I S (2-2) 

(with d the interlayer separation), or that it represents a short range interaction such that in 
momentum space it satisfies that V^/3(Q)Q^ vanishes at zero momentum. This includes the 
case of ultralocal potentials {i.e., with a range smaller or of the same order as the cyclotron 
length £), in which case we can set V^(0) = 0, or short range potentials with a range longer 
than I such clS cl Yukawa interaction. 

Following the same steps as in reference ^ , in the Appendix ^ we show that this system 
is equivalent to a system of interacting electrons coupled to an additional statistical vector 
potential a'^ (/i = 0, 1, 2) whose dynamics is governed by the Chern-Simons action 

Scs = T.^J d'x e^'^'^a^dA (2.3) 

a/3 

provided that the CS coupling constant satisfies 



27r(4siS2 - v?) 



2s2 —n 

y —72 2Si y 



(2.4) 



where Si, S2 and n are arbitrary integers. In eq (|2.3|) Xq.Xi and X2 represent the time and 
the space coordinates of the electrons respectively. In the equivalent theory the covariant 
derivative is given by 

D° = + i-^A, + lal (2.5) 

and it couples the fermions to the statistical gauge fields (a^), and to the external electro- 
magnetic field {A^. Notice that the theory has now a f/(l) ® U{1) gauge invariance. 

The Chern-Simons action implies a constraint for the particle density Jq (x) and the 
statistical flux = eijdia'j, given by 

jS{x)+k^(,B(,{x)=0 (2.6) 

(from now on we assume that repeated indices are contracted). 
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This relation states that the electrons in plane a coupled to statistical gauge fields with 
Chern-Simons coupling constant given by eq ( p.4| ), see a statistical fiux per particle of 27r2sQ, 
for the particles in their own plane , and a statistical flux per particle of 2?™ for the particles 
in the opposite plane. (Notice that in units in which e = c = % = 1, the flux quantum is equal 
to 27r). If the coefficient of the Chern-Simons term is chosen with the above prescription, all 
the physical amplitudes calculated in this theory are identical to the amplitudes calculated 
in the standard theory, in which the Chern-Simons field is absent. Of course, this is true 
provided that the dynamics of the statistical gauge fields is fully taken into account exactly. 

In this work we will take into account the dynamics of the Chern-Simons gauge fields 
in a semiclassical expansion, which is a sequence of well controlled approximations. In 
practice, we will only consider the leading and next-to-leading order in the semiclassical 
approximation. 

Using the constraint enforced by the Chern-Simons action, the interaction term of the 



action eq ( |2.1| ) becomes 

Sint = -\j d?z I dh' {k^'Bs{z) - -p^)V^p{\z- ^\){k''^B,{z') - -pp) (2.7) 

The quantum partition function for this problem is, at zero temperature 

Z[A,] = J VrV^Va''^ exp{zS{r, ^, A,)) (2.8) 

Since the action is quadratic in the fermions, they can be integrated out. The effective 
action (Ses) is given by (in units in which e = c = h = 1) 

5eff = -zEtr In {^^0 + + ^(^")'} + '5cs(a^ - i^) + 5:i^Ha^ - i^) (2.9) 

where 

S':i{a-, - i^) = -lJ dh J dh' {k'^\Bs{z) - Bs{z)) - p„) 

V^p{\z- ^\){kP\B,{z') - B,{z')) - -pp) (2.10) 

Here we have written the external electromagnetic field as a sum of two terms, one represent- 
ing the uniform magnetic field S, and a small fiuctuating term whose average vanishes 

7 



everywhere. The latter will be used to probe the electromagnetic response of the system. 
Notice that we have used the invariance of the measure Pa^ with respects to shifts, to move 
y4" out of the covariant derivatives and into the Chern-Simons and the interaction terms of 
the effective action (eq (p^)). 



A. Mean field approximation: allowed fluid states 

The path integral Z can be approximated by expanding its degrees of freedom in powers 
of the fluctuations, around stationary configurations of Sefj- This requirement yields the 
classical equations of motion. 

< jUz) >f = -f^ape,i[< £^{z) > -E^{z) >] 

-dteik J dh' n'^'V^siz, z') [-k'^' < B'' ~ W > {z') - p'] (2.11) 

Just as in the case of the single layer problem, these equations have many possible solutions, 
i.e., fluid states, Wigner crystals and non-uniform states with vortex-like configurations. 
We will only consider solutions with uniform particle density, i.e., the liquid phase solu- 
tion. This is the average field approximation (AFA), which can be regarded as a mean field 
approximation. At the mean field level the electrons in the layer a see a total flux B^g, 
equal to the external magnetic flux partially screened by the average Chern-Simons flux, 
i.e. B^ff = B -\- {B°') = B — (k~^)"^p^. It is easy to see that the uniform saddle-point state 
has a gap only if the effective field B^q is such that the fermions in layer a fill exactly an 
integer number pa of the effective Landau levels, i.e., those defined by B^f^. In other words, 
the AFA to this theory yields a state with an energy gap if the filling fractions of each layer 
satisfy 

1 

1^1 



±^ + 2si + 

Pi ' Ni 



1 

1/2 = — ^ (2.12) 



where A^i and N2 are the number of particles in layers 1 and 2 respectively. The sign in front 
of pi and P2 indicates if the effective field is parallel or antiparallel to the external magnetic 
field. 

Using the fact that z/ = z/i + z/2, and that the number of flux quanta enclosed by each 
plane is the same, the filling fractions can be written as follows 

n-i±^ + 2s2) 



n2-(±J^ + 2.i)(±^ + 2.2) 



n- (±^ + 2si) 

where pi, p2i -^i, ^2 and n are arbitrary integers. For the special case in which the two layers 
have the same occupancy, Ni = N2 and i^i = 1^2 = ^, the allowed fractions are 

where p is an arbitrary (positive or negative) integer. 

The effective magnetic field can be written in terms of the external magnetic field as 

^Sr = 5^ (2.15) 

For general values of these integers, the states whose filling fractions are given by eq ( ^.13[ ) 
have a gap at the mean field level of this approximation. This ensures that the perturbative 
expansion is meaningful. In other words, there is a small parameter, which is essentially the 
inverse of the mean field gap, for this perturbative expansion to be possible. If there is no 
gap for the excitations of the mean field ground state, the perturbative expansion breaks 
down. The breakdown is signalled by infrared divergencies at low temperatures. Such is 
the case for the compressible or "Fermi Liquid" states. We will see now that both types of 
states can occur for a given filling fraction (but not for all filling fractions!). 



It is clear from eq. ( |2.13| ) that in the bilayer systems, there are many possible choices 
of the numbers pi, p2, Si, S2 and n, i. e. , many different states, which have the same 
filling fractions. As a result, the phase diagram for bilayers is much richer than for spin 
polarized electrons in a single layer. Experiments on spin polarized 2DEG bilayers 
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at filling fractions v = 1,1/2 have shown a non-trivial phase diagram with at least two 
phases: a compressible phase in the regime in which the 2DEG's are well separated and an 
incompressible phase when they are closer by. However, as it was emphasized by Wen and 
Zee 1^, the incompressible states at these two filling fractions actually have quite different 
properties. We will see now that the phase diagram can be quite complex. 

Let us consider some cases of special interest, in particular with A^^i = N2- For example, 
at level one of the hierarchy, we choose pi = p2 = 1, 2si + 1 = mi, 2s2 + 1 = "012, and 
n = m, and we obtain the so called (mi,m2,m) states |l7|,|8|JT0(l whose filling fractions are 
u = "^^j^^^Sni" ■ -l-^ particular, for mi = m2 = 3 and m = 1, this is the (3, 3, 1) state, whose 
filling fraction is u = ^. For pi = P2 = 1, and 2si + 1 = 2si + 1 = n = m, we obtain the 



[m,m,m] states [[7|,p!T| whose filling fractions are p = 



We can also consider the limit = ^2 = P ^ 00. In this case, if Si = S2 = s, the 
state with filling fraction u = — will be obtained for all the values of n and s such that 
n + 2s = 2m. Since p —>■ 00, the effective field vanishes, and we find the analogous of the 
compressible states for the single layer problem discussed by Halperin et al |T^. All of these 
states are degenerate with the u = — states mentioned above. 

Clearly, it is always possible to construct a large family of, in principle, distinct states 
which have the same filling fraction. For instance, we can also choose pi = p2 = 1, n = —m 
and 2si + 1 = 2s2 + 1 = 3m. For the case Ni = N2, it is easy to see that all the states 
with the same value of n + 2s have the same filling fraction (at fixed pi = P2 = p)- For 
Ni 7^ N2, similar (but more complicated) families can be found. Let {Ni,N2) denote the 
largest common divisor {1. c. d. ) of A^^i and N2 and let i and k be two relatively prime 
integers (i. e. {^,k) = 1) such that A^i = i{Ni,N2) and N2 = k{Ni,N2). By direct 
inspection of the expressions for the filling fractions of eq. ( p.l2| ), we find that the following 
transformations leave the filling fractions Ui and 1^2 invariant: 

n ^ n — 2rk£ 
■Si — >■ si + r/c^ 
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S2^S2 + re (2.16) 

where r is an arbitrary integer. 

In simpler terms, new states are generated by transferring mtra-plane flux-particle at- 
tachment, determined by si and S2, to mier-plane flux attachment, determined by n. Since 
the inter-plane attached fluxes are negative in these states, the particles of one plane see 
the particles of the other plane as if they were holes. Thus, in these states, there is an 
effective attractive force between particles on different planes. Hence, in these states the 
wave functions, instead of having zeros of higher order (which represent repulsion), should 
have a larger weight when particles from different planes are closer to each other (in the 
sense of their coordinates projected on the xy plane). Asymptotically, these wave functions 
appear to have "poles" in the interlayer coordinates. Below we find that this is indeed the 
case. 

Which one of all of these states is realized for a given system at a fixed filling fraction 
should depend on the inter and intra layer interactions. At the mean field level we find that 
these states are all degenerate. In appendix ^ we give an estimate of the contribution of the 
gaussian corrections to the ground state energy. Our results show that this degeneracy is 
lifted and that, of all the incompressible states at filling fraction 1/m, the simpler (m, m, m) 
states have lower energy. Nevertheless, the situation should be in principle more complex. 
For instance, one can imagine an interlayer interaction which is very weak at short in- 
plane distances but stronger at intermediate separations ("hollow core"). In such cases, 
the off-plane attachment mechanism might prevail. In addition, for truly large interlayer 
separations, the system should prefer to have (for v = 1/m) a state with a "Fermi Liquid" 
state on each layer. A detailed calculation of the contribution of fluctuations to the ground 
state energy is necessary to establish which one of these states actually occurs for a given 
form of the interaction and interlayer separation. We will return to this problem elsewhere. 

As it stands, this theory has no a priori way of limiting the number of possible states. 
Clearly, many of the states in these families must be redundant since, for a fixed number 
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of particles and fluxes, the dimension of the Hilbert space is flxed. Hence, either by virtue 
of a symmetry the states in these famihes faU into equivalency classes, or there should be a 
natural way to limit the number of states. Fluctuations must also play an important role 
here. Already at the gaussian level the fluctuations do lift these degeneracies. However, at 
the present time, it is unclear how will the fluctuations manage to reduce the size of the 
Hilbert space. 

A more natural resolution of these issues would be that a physical interpretation of these 
new families of states could be found in terms of the physical degrees of freedom of the 
electrons. It is natural to speculate that the degeneracy found here may become natural 
in the spin picture of the bilayer system. Thus, if we think the electrons in the upper 
layer as having spin up while those in the lower layer have spin down, the states can now 
be classified in terms of the total number of electrons (charge), and of the difi'erence of 
occupancy, which can be viewed as the z-projection of the "spin". All of the states 
that we have discussed here have the same occupancy in both planes and, hence, they 
have Sz = 0. In the limit of small separations, the system should have an effective SU (2) 
symmetry even if the electrons are fully polarized. In that limit, there is an additional 
operator which commutes with the Hamiltonian, S^, the "total spin". Even though for an 
arbitrary interlayer interaction the SU{2) symmetry is broken down to the SU{2) 

states can still be used to enumerate the states (even though they are no longer good 
quantum numbers). Thus, it is natural to think that the multiple solutions may be linear 
combinations of the = ^ states with different total "spin" S < 2N. For the counting of 
states to be the same, it is necessary that the families of states found above either terminate 
(and be finite) or become split into equivalence classes. In this system, the 5'C/(2)-invariant 
limit is achieved either when the two layers physically coincide or if the electrons are non- 
interacting. However, the two-component Chern-Simons gauge theory that we use here does 
not exhibit the SU(2) symmetry explicitly even in the limit of zero interlayer separation. 
In this theory, the SU{2) symmetry is a dynamical symmetry. In other terms, given that 
this is an exact transcription of the two-component 2DEG, the exact states should form 
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SU{2) multiplets. However, while this is a property of the exact theory, it is by no means 
guaranteed that a semiclassical approach would be able to recover this exact property. Such 
dynamical mechanisms are known to work in the abelian bosonization of SU{2) invariant 



1 + 1-dimensional Fermi systems ||13|| . 



Finally, since the bilayer system can also be used to describe spin unpolarized electrons, 
the allowed bilayer states should include partially polarized and unpolarized (singlet) spin 
states. One problem in the way of making this connection is the fact that this formalism 
breaks down if the Chern-Simons coupling matrix is singular. This happens for n = ±2s. 
Thus, the spin singlet state (3, 3, 2), which has filling fraction u = 2/5, cannot be described 
within this abelian Chern-Simons approach. In the section |VI1| we present a generalization 



of the theory of Balatsky and Fradkin which describes the SU (2) symmetric cases. 



B. Semiclassical approximation (RPA) 

Now we consider the gaussian (or semiclassical) fluctuations of the statistical vector 
potential around the mean-field state. The gaussian corrections must alter the qualitative 
properties of the state described by the AFA, which violates explicitly Galilean invariance 
(more generally, magnetic invariance) which, for translationally invariant systems, must 
remain unbroken and unchanged. Thus the center of mass of the system must execute a 
cyclotron-like motion at, exactly, the cyclotron frequency of non interacting electrons in the 



full external magnetic field, as demanded by Kohn's theorem [|T^]. A naive application of 
the AFA would suggest that the cyclotron frequency is renormalized downwards since the 
effective field seen by the composite fermions is smaller than the external field B. Hence, 
the magnetic algebra may appear to have changed. In the same way as it happens for the 
single-layer systems |TB[ , here the gaussian fluctuations yield the correct cyclotron frequency 



and, thus, restore the correct magnetic algebra. 

At the gaussian level, the effective action for a° is 
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-^1 d'x I d'y {Bs{x)-Bs{y))K'^'V^p{\x-y\)K^^{B,{y)-B,{y)) 

+ ^ / d'x e^^\dl - ~Al)dAd{ - ~Ai) (2.17) 

The tensor n^^(x, y) = Yl^{x^ y)Sa/3, where U'^^^x, y) is the polarization tensor of the equiv- 
alent fermion problem at the mean field level. It is obtained by expanding the fermion 
determinant up to quadratic order in the statistical gauge field. This tensor was calculated 
in reference 0|. The subindex a indicates that the effective field which appears in the 
expressions for li^^ in reference is i?cfr = 

After integrating out the gaussian fluctuations of the statistical gauge fields a^, we obtain 
the effective action for the electromagnetic fluctuations (A^), Sl^^A"^) 

SmK) = ^ / / y)At{y) (2.18) 

Here K'^^ is the electromagnetic polarization tensor. It measures the linear response of the 
system to a weak electromagnetic perturbation. 

We will use this effective action to calculate the full electromagnetic response functions 
at the gaussian level. Since this calculation is based on a one loop effective action for the 
fermions {i.e. a sum of fermion bubble diagrams), this approximation amounts to a random 
phase correction to the average fleld approximation. 

The components of the electromagnetic polarization tensor can be written in momentum 
space as follows 

K^f = uQ.Kficu, Q) + le.kQkK ficu, Q) 
K"^ = ujQjKo'^iuj, Q) - iejkQkKi'^ii^, Q) 

Ktf = u^6,jKfiu, Q) - ie,,uKf{u, Q) + {Q% - Q,Q,)Kfiu, Q) (2.19) 
where K"^{uj,Q) {i = 0, 1,2) can be written as a power series expansion in powers of 

cff 

and have poles for the values of u that coincide with the collective modes of the system. 
The electromagnetic response functions determined by K^^j^ have the following properties: 
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i) As in the single-layer case, the collective excitations of this systems are determined by 
the poles of the density correlation function, Kqq{uj,Q). 

ii) The leading order term in of the K^q component of the polarization tensor saturates 
the /-sum rules. These sum rules correspond to the conservation of the particle number in 
each layer separately. This result is essential in order to show that the absolute value squared 
of the ground state wave functions of all the (incompressible) liquid states have the form 
described in the introduction at very long distances and in the thermodynamic limit. 

iii) The gaussian fluctuations of the statistical gauge field are responsible for the FQHE. 
In particular, the gaussian corrections yield the exact value for the Hall conductance. 

In the next sections we will discuss these properties in detail. 

III. SPECTRUM OF COLLECTIVE EXCITATIONS 

In this section we derive the spectrum of collective excitations for two different states, 
the (m, m, n) and the (m, m, m) states. We use the same approach as we did for the single- 
layer case [^], i.e., we study the poles of the density correlation function to determine the 
collective modes of the system. 

A. Collective excitations for (m, m, n) states 

For simplicity, we have studied the collective modes for the state (3, 3, 1). All the other 
states can be studied by straightforward application of the same methods. 

In this case the total filling fraction is z/ = i, being ui = U2 = \ . The effective cyclotron 
frequencies and magnetic fields are ujI^ = cj^g = uj = ^ and Bl^ = B^^ = B = ^. 

We find that there is a family of collective modes whose zero-momentum gap is kO, where 
k is an integer number different from 1. At mean field level, there are two modes for each 
integer multiple of Co. After including the gaussian fluctuations we find that there are no 
modes with a zero momentum gap at u. One of them has been "pushed up" to the cyclotron 
frequency and the other up to 2uj (at Q = 0). Therefore, at these multiples of Co there are 
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three degenerated modes for Q — 0. For Q 7^ 0, the degeneracy is hfted and these three 
modes have different dispersion curves. 

At 2uj we find that there are two modes with residue Q^, and one with residue Q'^. The 
former are 



2-Befr 



-2 

J2 o; a± 



(3.1) 



where 



3M, 

a± = -t;—{vu - V12) ± 

ZTT 



(3.2) 



Here v^/j are the zeroth order coefficients of the Fourier transform of the interparti- 

2 

cle pair potential for short range interactions. For Coulomb interactions vu — ^ and 
V12 = ^e-l^l'' Ri ^ if \Q\d < 1, therefore a± = ±4 in this limit. 
The residues in Kqq corresponding to these poles are 



Stt 



1 -1 
-1 1 



(3.3) 



It is clear from the form of these residues, that these excitations are out of phase modes, 
because they only couple the "out of phase" density (i.e., they couple p~ with itself). 
The other mode at 2uj is 



c^o((3)= (2cu)'-6( 



2 al ^ ^ -2 



2B, 



eff 



(3.4) 



and its residue is proportional to 



li i\ 

1 1 



Thus, this is an in phase mode since it only 



couples p+ with itself. 

The two modes with zero momentum frequency Sou are given by 



(3.5) 



/ 



with residue proportional to 



1 -1 
-1 1 



(this is an out of phase mode); and 
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2B. 



eflf 



(3.6) 



with residue proportional to 



(i i\ 



(this is an in phase mode) . 



There are three modes whose zero momentum frequency is the cyclotron frequency. 



with residue 



u{Q) = u;, + (2 + ^ 



TT 



V 



Res{KoQ,oj{Q)) ^ -Q ujc— 



Stt 



(1 i\ 



1 1 



(3.7) 



(3.8) 



2 iw Q 



2-Befr 



^2 _2 ^ ^ 3M(wn^^^_i 



27r 



with residue proportional to 



^ 1^ 



1 1 



(3.9) 



-0^ + 2(^)^0;^ 



with residue proportional to 



( 1 -i\ 



-1 



(3.10) 



. The first two modes at cUc are in phase and 



the last one is out of phase. 

In summary, we find a family of collective modes with dispersion relations whose zero- 
momentum gap is kij, where k is an integer number different from 1. When k = 4, i.e. the 
zero-momentum gap is the cyclotron frequency, there is a splitting in the dispersion relation 
for finite wavevector. However, only one of these three modes has residue proportional 
to in the density correlation function. One expects that the other modes will become 
damped due to non-quadratic interactions among the collective modes. On the other hand, 
the collective mode with lowest energy which has A; = 2, is stable (at least for reasonably 
small wavevectors) . 
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The validity of the spectrum described in this section is hmited by the fact that we 
have not considered the physics at arbitrary wavevectors, and the (expected) effects of non- 
gaussian corrections. At the gaussian (RPA) level and for small momentum, we found a 
family of collective modes which are infinitely long lived, {i.e., the response functions have 
delta- function sharp poles at their location). These modes represent charge- neutral bound 
states. For Q sufficiently large, the energy of the collective mode can become equal to the 
energy necessary to create the lowest available two-particle state: a quasiparticle-quasihole 
pair. At this point, the collective mode should become damped. Non-gaussian corrections 
to the RPA are also expected to give a finite width to (presumably) all the collective modes 
but the lowest one. Since the modes with zero momentum gap at kCu, k > 3, are not the 
collective modes with lowest energy, it is possible that at finite wavevectors they may also 
decay into the collective mode with lowest energy ( the mode with k = 2, which has a gap 
at cu). 

B. Collective excitations for {m,m,m) states 

Here we present the spectrum of collective excitations for the so called (m, m, m) states. 
In this case the total fiUing fraction is u — being ui — U2 — . The effective cyclotron 
frequencies and magnetic fields are cj^jf = cj^g = Co — and B^^^ = B'^g = Bes = 

We find again a family of collective modes whose zero-momentum gap is kuj, where k is 
an integer number different from 1. At mean field level, there are two modes for each integer 
multiple of a). After including the gaussian fiuctuations we find that there are no modes with 
a zero momentum gap at uj. One of them has been "pushed up" to the cyclotron frequency. 

— * — * 

Therefore, at Uc there are three degenerate modes for Q = 0. For Q ^ the degeneracy 
is lifted and these three modes have different dispersion curves. The other mode at u has 
been "pulled down" to zero frequency at Q = 0, i.e., it has become a gapless mode. 

We will distinguish between the cases m — 1 and m ^ 1. 
a) Case m — 1 
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The gapless mode is 



where 



u;{Q)^Vs\Q\ 



M , 



(3.11) 



UJcV 



2M 



(3.12) 



where are the zeroth order coefficient of the Fourier transform of the interparticle pair 

— * 2 — * 

potential for short range interactions. For Coulomb interactions v\\{ff) — ^ and vviiff) — 



-\Q\d 



if <^ 1. Therefore, {v\\ — Vy^ — for Coulomb interactions (in the limit 



\Q\d^\). 

The residue in K^^ corresponding to this pole is 



Stt 



V 



-1 1 



(3.13) 



Therefore, this is an out of phase mode. 

At ujc — 2uj we find that there are two (m phase) modes with residue 



^±(g) 



^1 + 



- - 2 

1 2 a; a± 



(3.14) 



where 



M / M 



(3.15) 



For Coulomb interactions Vn + V12 — '2^ li \Q\d <^ 1, therefore this term is higher order in 
Q and it should be neglected, i.e., a± — 4. 

The residues in Kqq corresponding to these poles are 



Res{Koo,(jJ±{Q)) = -Q ojc— 

OTT 



1 1 



(3.16) 



therefore these are in phase modes. The other mode at uJc is 

2 



-2 
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)a; 



(3.17) 



and its residue is proportional to Q"^ 



1 -l\ 



V 



-1 1 



[out of phase mode). 



The modes with zero momentum frequency uo = kuo for A; > 3 coincide with the expres- 
sions given below for the case m 7^ 1, in eq ( p.22|) and eq ( ^.231) . 

b) Case m 7^ 1 

The gap less mode is an out of phase mode with the same form as for m = 1 (eq ( |3.11| )) 
and with the same residue (eq ( |3.13| )). 

At ojc = 2mLJ we find that there is one in phase mode with residue 



+ — (t^ii + t^i2) 2m(-|-) 



m — 1 271 
The residue in Kqq corresponding to this pole is 



V 



57r 



1 A 



The other modes at lOc are one in phase mode 

4m2(2m-l)2(2m-2) 



\2m-2 -2 



(2m - 2)! ((2m - 1) + (m - l)§{vu + ^^12) 



'2R 



00 



cff 



with residue proportional to Q'^^"^ ^\ and one out of phase mode 



a;(Q) 



I \2m-l -2 



(2m-2)!'25eff' 
with residue proportional to Q^"^. 

The two modes with zero momentum frequency kuj {k ^ 1, 2m) are given by 

2 2 , 



uj{Q) 



kuj 



(A;-2)!'25eff 



and 



kuj'~ 



2k{2m-l) g2 ^ 



Their residues are proportional to Q 



{k -2)\{2m- k)^2B^s 

( 1 -i\ 

1 1 



2k 



\ 



, and to Q 



2k 



I 



(l l\ 

1 1 



(3.18) 



(3.19) 



(3.20) 



(3.21) 



(3.22) 



(3.23) 



respectively. 
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In summary, for the (m, m, m) states, in addition to a family of collective modes with 
dispersion relations whose zero-momentum gap is ku {k ^ 1), we find that there is a gap less 
mode. The gapless mode is related to the relative fluctuations of the electronic density in 
each layer, i.e., to the fluctuations of — ■ 

All the considerations about the validity of this spectrum of collective excitations beyond 
the semiclassical approximation that we discussed in the previous section are of course valid 
in this case. 



IV. GROUND STATE WAVE FUNCTION 

We need to show first that the long wavelength form of Kqq, found at this semiclassical 
level, saturates the /-sum rule. This result implies that the non-gaussian corrections do not 
contribute at very small momentum. We will use this result to show that the absolute value 
squared of the ground state wave function of this state has the Halperin |0] form at very 
long distances, in the thermodynamic limit. 

The /-sum rule can be derived as follows. The retarded density and current correlation 
functions of this theory are, by definition 

Df,{x,y) = -ze{x,-y,)x<G\[r^{x),jP^y)]\G> (4.1) 

where (/x = 0, 1, 2) are the conserved currents of the theory, and > is the ground state 
of the system. Using this definition and the commutation relations between the currents, 
one can derive the /-sum rule for the retarded density correlation functions Dqq. In units 
in which e = c = % = 1, it states that 

11 ^^'^^oo'M)='^Q'6-^ (4.2) 

This equation implies the conservation of the current in each layer separately. It is easy 
to show that, in the basis of the total and relative currents, i.e., in the p± = p^^p^ basis, 
eq ( [4. 21 ) states the conservation of p+ and p_ independently. 
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On the other hand, it can be shown (see for instance reference [|T9[) that the polarization 
tensor K'^^ and the density and current correlation functions V^^^ satisfy the following 
identity 



K?:^{x,y) = -Df^{x,y)+< 



> 



(4.3) 



Thus, eq ( [4 .21) is also valid if we replace D'^£ by Kqq and we change the sign in the r.h.s. 



of the equation. 



A. Ground state wave function for (m, m, n) states 



We have found that for an (m, m, n) state, the leading order term in of the zero- zero 
component of the electromagnetic response is given by 



P 



AM u;2 



P 



uf. + It 



^ 1^ 



1 1 



AM uj'^ - ((m - n)ujY + ie 



I 



\ 



-1 1 



(4.4) 



Notice that the poles whose residues in K^q are proportional to Q^, have zero momentum 
frequency given by (m + n)u! = ujc and (m — n)uj. In particular, if m = 3 and n = 1, i.e., 
for the (3, 3, 1) state, the expression in eq (|4.4| ) coincides with the result given by eq ( p.3|) 



and (|3.8|) , provided that 



Tf , and that (1 + ^) ^ + (1 



16 \-l 



1, with a±. defined 



by eq (|3.2|) . In the gaussian approximation to the Chern-Simons Landau-Ginzburg theory 
for the double layer systems |17|JTT[], these two modes have also residue proportional to 
in the density correlation function. But in that approach, each of these modes separately 
saturates the /-sum rule. 

The correlation functions that we derive from the path integral formalism are time- 
ordered. Therefore, if we use the relation between time-ordered and retarded Green's func- 
tions, and eq ( |4.3D and (|4!^ ), we see that the leading order term of K'^^ saturates the /-sum 



22 



rule, eq ( |4.2|) , already at the semiclassical level of our approach. Notice that for this state 

p' = p' = i 

Thus, the fermionic Chern-Simons approach gives the correct leading order form for the 
density correlation function, in the sense that it is consistent with the /-sum rule, at the 
semiclassical level of the approximation. 

It is important to remark that the coefficient of the leading order term of Kqq can not 
be renormalized by higher order terms in the gradient expansion, nor in the semiclassical 
expansion. In the case of the gradient expansion, it is clear that higher order terms have 
higher order powers of Q^, and then, do not modify the leading order term. In the case 
of the corrections to Kqq originating in higher order terms in the semiclassical expansion, 
they also come with higher order powers of Q^. The reason of that is essentially the gauge 
invariance of the system. This implies that the higher order correlation functions must be 
transverse in real space, or equivalently they have higher order powers of in momentum 
space. Being higher order terms in the expansion they can not change the leading order 
term. 



We will now follow the method used in reference |2y] to write the ground state wave 
function in the density representation. We begin by recalling that the absolute value squared 
of the ground state wave function in the density representation |\E'o[p]P is given by |2l[ 



|^o[pi,P2]|' = /^^^o e-'/'^'"^o(-)p"(-) ^ lim^ (0|Te*/'=''^' ^?(-) jo"W|o) (4.5) 

where Jq (^) = P^i^)- The operators in this expression are Heisenberg operators of the 
system in the absence of sources. The vacuum expectation value in the integrand of eq ( |4.5|) 
can be calculated from the generating functional of density correlation functions, Z[A'^] . 
Clearly, we have 

\MpuP2]\' = I -DA-o e-/'^^^^o(^Tp'^(^T lini Z[A-,,A- = 0]. (4.6) 

Eq ( ^^ ) tells us that | \E'o [pi , P2] P is determined by the generating functional of equal-time 
density correlation functions. 
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The generating functional that appears in that expression is given by 

hm 2:[A^,A" = 0] = / Vi)*Vi)Val e^^^^*''^'""'^") (4.7) 

A'^{x)^A'^{x)5(xo) J 

The path integral on the r.h.s. of eq ( |4.7| ) can be written in terms of the effective action 
S(.ff{A^) for the external electromagnetic field. We have seen that, in the thermodynamic 
limit, and for weak fields, the effective action admits the expansion given by eq ( |2.18| ). Since 



we need only the density correlation functions, it suffices to know the zero-zero component 
of K'r^^. In momentum space, and in the small limit, K^q is given by eq ([4.41). We can 



see that the dominant term in Kqq is of order \/B. Higher order terms in the gradient 
expansion will contribute with higher powers of 1/B. The same observation applies for all 
the corrections to Kqq originating in higher order terms in the semiclassical expansion. Here 
the thermodynamic limit is crucial since we are only taking into account fluctuations with 
wavelengths short compared with the linear size of the system. The higher order terms, 
which vanish like powers of /B^ can only be neglected for an infinite system. 
Using eqs (p.l8|) and ([4. 41), eq ([4.7|) becomes 



lim Z\An A" = 01 = e^/'''^'^''^^o(^)(ii'^-o-»TO-^M (^•'S/))^o(y) (4 8) 

or, by Fourier transforming the exponent, we get 

2 °° 

lim = 0] = e — (4.9) 

A«(x)^A-{a;)5(xo) 



The terms dropped in the exponent of eq (|4.^ ) and eq ( |4.9| ) represent equal-time density 
correlation functions with more than two densities. These terms give rise to three-body 
corrections ( and higher) to the wave function and modify the Jastrow form. The kernels 
of these non-linear contributions are, by gauge invariance, required to be transverse. Thus, 
in momentum space, the residues of their poles have higher powers in 

Q2 than K^^{uj,Q). 

Since, by dimensional analysis, each power of has to come with a factor of these 
terms which are not bilinear in the densities are subleading contributions in the limit -B — > 00. 
At the level of the gaussian ( or semiclassical) approximation, these kernels are equal to zero. 
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All of these considerations hold provided that the Fourier transform of the pair potential 
satisfies Q^V{Q) as ^ 0. 

Replacing the expression for Kqq{uj, Q) given by eq ( [4 .41 ) into eq ( [4.9[ ) and integrating 
out Aq, we obtain the following form for the absolute value squared of the wave function 

Ni N2 Ni N2 

= n i^i-^jf" n i^i - %f n n - %f " 

i<j=l i<j=l i=l j=l 

n Ni N2 

exp{-^(Ei^^r+Eb"'^r)} (4-10) 
^ i=i i=i 

where the coordinates Xi are in plane 1, jji are in plane 2, and Ni = N2 = y this state. 

To get this result we have used that the eigenvalues of the local density operator, in 
a Hilbert space with particles, are p°'{x) = Yli^=i^i^ ~ ^i) ~ P°- A similar result was 
obtained recently by Schmeltzer and Birman |22| who used a different approach. 

Notice that the wave function of eq ( [4.10| ) is the absolute value squared of the Halperin 
wave function Numerical calculations have established ||2^ that this wave function 
accurately describes the ground state wave function for the (3, 3, 1) state for d = 1.5ic, 
where ic is the cyclotron length. We have shown that eq (|4.1(]| ) gives the exact form of the 
ground state wave function at long distances and in the thermodynamic limit. Since the 
leading order term of Kqq saturates the /-sum rule, higher order corrections in the expansion 
cannot modify this result. 

In the same way as for the single layer problem, we have shown that, in the thermody- 
namic limit, the exact asymptotic properties of the wave function, when its arguments are 
separated by distances long compared with the cyclotron length (but short compared with 
the linear size of the system), are completely determined by the long distance behavior of 
the equal-time density-density correlation function [i.e., the structure factor). 



B. Ground state wave function for (m, m, m) states 



For these states, the leading order term in of the zero-zero component of the electro- 
magnetic response is, according to eq (|3.13|) and ( p.l6|) or ( p.l9|) 
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where we have used that ■§ = and that [(1 + ^) ""^ + (1 + ^) ""^l = 1 with a± defined 



by eq§J^. 



Following the same steps as in Sec |IV A| we can prove that the leading order term of 
Kqq saturates the /-sum rule, eq (|4.2|) . All the remarks about the exactness of this result 
are also valid in this case. In this case too, the Chern-Simons Landau-Ginzburg approach 
violates the /-sum rule at the gaussian level |17|JTT[]. 



Substituting eq ( ^4.11|) into the expression for the absolute value squared of the ground 
state wave function (eq ( [4.6| ) and ( |4.9| )) we obtain 

A^i N2 Ni N2 

i^(fi,...,xivi,yi,...,f7vji' = n n nn i^i-^ii"" 

i<j = l i<j = l 1=1 j = l 

TJ Ni N2 

i=l i=l 

\^ij^i\Xt Xj\ ij=i\yt Vji i=ij=iFi Uji J 

(4.12) 

where the coordinates Xj are in plane 1, yi are in plane 2, and Ni = N2 = y ^^r this state. 

Notice that this ground state wave function is not exactly the same as the (m, m, m) 
Halperin wave function, to which the true ground state approaches as d 0. There is 
an extra factor which comes from the fact that there is a gapless mode in the spectrum 
of collective excitations. This contribution is analogous to the phonon contribution to the 
wave function of superfluid He4, and just as in that problem, it is essential to obtain the 
correct properties for the spatial correlations of the ground state. This contribution is very 
small at long distances compared to the cyclotron radius, and it can not be written only in 
terms of coordinates in the lowest Landau level. 
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Based on the same arguments that we discussed in Sec [IV A| , we can argue here that 
eq ( |4.12| ) is an exact resuh for the asymptotic form of the ground state wave function , at 
long distances and in the thermodynamic hmit. 



V. HALL CONDUCTANCE 

We show now that, already within our approximation, this state does exhibit the Frac- 
tional Hall Effect. In order to do so, we will calculate the Hall conductance of the whole 
system. 

Since we are only interested in the leading long-distance behavior, it is sufficient to keep 
only with those terms in the action of eq ( ^.17] ) which have the smallest number of derivatives, 
or in momentum space, the smallest number of powers of Q. Therefore, the leading long 
distance behavior {i.e., small momentum) of the effective action for the fluctuations of the 
Chern-Simons gauge fields and of the electromagnetic field is governed by the Chern-Simons 
term. In this limit eq ( p.l7|) turns out to be 



where = u and = —Qi according with the convention that we have used in reference 
@, and R^^ = + 

The next step is to integrate the statistical gauge fields to obtain the effective action for 
the electromagnetic field. In particular, we will need to compute the inverse of the matrix 
^a/3_ This inverse only exists if A = [(^ + 2si)(^ + 2^2) — n^] 7^ 0. Therefore we must 
consider two cases, the one in which A 7^ and the one in which A = 0. 
Case A ^ 

Upon integrating over the statistical gauge fields in eq the effective action for the 

electromagnetic field results 

^eT(^^) - - {---'V'-''] J e,^xQ' Al (5.2) 
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where the coefficient = k'^'^[5'^^ — {k ^)^^k^^] is given by 



a/3 
^eff 



27r[(J^ + 2si)(^ + 2.2)-n2 



V 



-n 



-n 



(5.3) 



In particular, if we consider the case in which both layers are coupled to the same electro- 
magnetic field, then Aj^ = Aj^ = A^, and the coefficient in the effective action results 

1 2n-(J^ + 2.i)-(i7 + 2s2) 



E 



.a/3 
"eff 



•-P2 



V 

2^ 



(5.4) 



27r n2 - (^ + 2si)(^ + 2s2) 

The electromagnetic current induced in the system is obtained by differentiating the 
effective action S'efT(A^) with respect to the electromagnetic vector potential. The current 
is = ^e^ivxP"^ ■ Thus, if a weak external electric field Ej is applied, the induced current 
is Jfc = KcS^ikEi- Therefore the coefficient K^g is the actual Hall conductance of the system. 



xy 



V 

Kpflp = 

27r 



(5.5) 



which is a fractional multiple of ^ (in units in which e = h = 1). Thus, the uniform states 
exhibit a Fractional Quantum Hall effect with the correct value of the Hall conductance. 
Case A = 

It can be shown that when is not invertible, i.e., it has a zero eigenvalue, the cor- 
responding linear combination of the gauge fields become massless. In other words, the in 
phase gauge field = + has a finite gap which couples to the electromagnetic field 
Afj_, while the out of phase gauge field ci~ = a^^ — is gapless. 

We will study in particular the case of the (m, m, m) states which satisfy the condition 
A = 0. For these states 



—m 



K 



27r(l -2m) 



(5.6) 



which is clearly non invertible. 

We can write the effective action for the fluctuations of the Chern-Simons gauge flelds 
(eq ( |5.1] )) in the basis deflned by = ± . 
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(27r)3 


^ 27r(l 


— 2m 


d'^Qduj 


1 
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(27r)3 


47r(l - 


2m) ^ 


d^Qduj 


1 




(27r)3 


47r(l - 


2m) *^ 


d^Qdu 


1 





The gauge field d~ appears as a Lagrange multiplier in this action. The integration over 
it states that the current e^^xQ^ Ay vanishes. This is trivially valid if the electromagnetic 
field is the same for both layers, because A~ = 0. 

The integration over a+ gives the effective action for the field A^. If we consider the 
case in which A], = Af, = An, the result is 



d^Qduo ~ 



Following the same steps as in the previous case, the Hall conductance results a^y = = ^ 
which is the correct value for the (m, m, m) states. 



VI. QUANTUM NUMBERS OF THE QUASIPARTICLES 

In this section we evaluate the charge and statistics of the quasiparticles. We will closely 



follow the methods and notation of reference 19 



We need to identify the operator which creates the quasiparticles within the framework 
of the Chern-Simons theory. Let us consider the gauge invariant operator which creates an 
excitation at the point at t = in layer a, and destroys it at the point y at time t = T in 
the same layer, and which behaves as a quasihole 

Here r{x,y) is a path in space-time going from {T,y) to {0,x). In this expression we are 
defining the quasihole operator in the layer a, therefore the index a is fixed, i.e., there is no 
sum over a assumed. The operator in eq (|6.1| ) is invariant under gauge transformations of 
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the statistical gauge field, but depends on the choice of the path F. In this expression, the 
fluctuations of the electromagnetic field have been switched off. The system only feels 
the uniform magnetic field determined by A^, and the statistical gauge fields. 
Our goal is to evaluate the Green function Gr{x,y) defined by 

Gr{x,y) = {GS\T [C(x)e(-*^r,.,„K+^^)'^-^) (6.2) 

where T is the time ordering operator. 

We calculate this Green function in the path integral formalism, where it is given by 
an average over the histories of the fermionic and statistical gauge fields, weighted with the 
amplitude exp{iS{ip*,ip,ci'^,A^)), where the action is defined by eq ( p.8|) . After integrating 
out the fermionic fields, the Green function can be written, up to a normalization factor, as 

Gr{x, y)= I Val G«„(x, y) e-^i'r(.,.)K+^M)<i-M ^g g) 



where S^s is given by eq (|2.9|) . The function Gaa{,x, y) is the one particle Green function for a 
problem of fermions in a field determined by the statistical gauge field (a^) plus the external 
magnetic field {A^), and at finite particle density determined by the chemical potential 
(/io). The one particle Green function can be written in terms of a Feynman path integral 
as follows 



G,„(x, y) = -te^^^ I Vz[t] e'^^'^'^^ (6.4) 
with the boundary conditions 

lim2;(t) = X 

\m^z{t)=y (6.5) 

The weight e^*^^-* serves to fix the number of particles. Since the mean field solution has 
effective Landau Levels filled, the chemical potential has to be set to lie between the levels 
•pa and + 1. The path integral in eq ( |6.4D is a sum over all the paths F which go from 
X to y in time T. The action is the standard one for non-relativistic particles coupled to a 
gauge field 
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Since there is an energy gap in this problem, in the long-distance, long-time limit, the path 
integral is dominated by paths close to the solution of the classical equations of motion. 
Therefore, the dominant trajectories are smooth. Thus, it should be a good approximation 
to pull the integral over the trajectories of the particles z{t) outside of the functional integral 
over the statistical gauge fields. The integral over the trajectories will be done at a later 
stage. We can write the Green function as 

Grix, y)^l Vz[t] e'lo f (f )^ J Va^ e"^ ^^^2+^.)'^-. ^^s., (g 7) 

In this expression, the set of closed curves 7 represents paths which are the oriented sum 
of the path F from y to x (which is fixed), and all the possible paths 7 from x to y (which 
correspond to the histories of the particles). 

In the semiclassical approximation, the path integral over the statistical gauge fields 
is replaced by an expansion around the solutions of the classical equations of motion. In 
this approximation, the particle only feels the electromagnetic field screened by the average 
of the statistical gauge fields. In other words, the effective field felt by the particles is 
B^Q = B — (k~^)"^p^. It is clear from eq (|6.7| ) that, for each closed trajectory 7 there is a 
constant factor which can be factored out from the functional integral, and that corresponds 
to an Aharonov-Bohm phase factor for a particle moving in the field 5"^, not in the external 
field B. In fact, the exponent of the Aharonov-Bohm phase factor is equal to B2gA±{'y), 
where ^±(7) is the (spatial) cross sectional area bounded by the path 7. Defining the 



effective charge as = and using eq (|2.15|) , we find that 



^eff = ^ (6.8) 

In particular, for the (m, m, n) states, the effective charge in both layers is the same and is 
given by qes = (^n+m) • This result coincides with the one in reference |p5|| . 

The fractional statistics can be studied by considering the two particle Green function. 
The generalization of the above formalism for this case is straightforward. The only difference 

31 



is that, for the two particle case, there will be two sets of trajectories, one for each particle. 
We will discuss first the case in which both particles are in the same layer, let say, the layer 
a. The two particle Green function Gp ^(xi, 0:2, yi, 2/2) is defined by 

G'f^^(Xi,X2,l/l,2/2) = 

(6.9) 

Following the same steps as for the one particle case, we can write this two particle Green 
function in terms of a path integral over the statistical gauge fields and over the trajectories 
of the particles. In this case, the Grassman integral automatically antisymmetrizes the two 
particle Green function, and it comes as a sum of direct and exchange processes with the 
gauge fields as a fixed background. The two particle Green function turns out to be 

where the path 7^ corresponds to direct processes (where particle 1 is destroyed at X\ and 
created at and particle 2 is destroyed at and created at 1/2), and the path 7e corresponds 
to exchange processes (where particle 1 is destroyed at X\ and created at yi, and particle 2 
is destroyed at x^ and created at y\). Note that there is a relative sign between these two 
processes. 

In the low energy limit, the dominant paths are very long and wide. Therefore, to 
compute the integral over the statistical gauge fields it will be sufficient to consider the 
effective action in the infrared limit. This effective action only contains the Chern-Simons 
term, and it is given by eq ( |5.1| ) but taking = 0. For both, direct and exchange processes, 
we need to calculate averages of the form (exp[i dx'^a'^])cs, where the subindex CS indicates 
that we only keep the Chern-Simons term in the effective action. If the coefficient k"'^ in the 
effective action is invertible, we can follow the steps described in Appendix |^ to calculate 
these averages. The result is 
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Recall that the index a is fixed and indicates the layer to which the quasiparticle belongs. 
The current j° is a three-vector of unit length tangent to the world lines and takes a non-zero 
value only on the world lines of the particles (7). The integral §^ da nJJ j° counts the number 
of times the current pierces the surface a, therefore, it is equal to the linking number 
of the curve 7, z/^. The configuration of paths can be classified according to their linking 
number. The weights of the configurations with different linking numbers have different 
phase factors. Also, configurations of paths from direct and exchange processes also have 
different linking number. While the phase factors themselves depend on the trajectories, and 
thus on the arbitrarily chosen paths for the two particles, the relative phase only depends 
on the topological properties of the configurations of paths and it is determined entirely by 
the relative linking number Au^. In particular we compare two paths which form a linked 
knot with two paths which do not. In this case Az/^ = 1. Therefore, the difference of phase 
between the direct and the exchange terms for these kind of paths in the two particle Green 
function is 

(j) J ^ ) 

S^ = n(l + . — '^^^r" ) (6.12) 

where if a = 1, then (3 = 2 and viceversa. In particular, for the {m,m,n) states, the 
statistics is 5 = — (^2'!^ ^2) , independent of the layer. This result coincides with the one in 
reference ||25|| . 



We define now the relative statistics of two particles in different layers, as the relative 
phase factor that we obtain if we compare the following processes 

{GS\T[iljl{xi) ^2*(x2)e^"'^r(.,..i)K+^-)'='---*/r(.2,.2)('^"+^^)''"'^) ^i{yi)My2)]\GS) (6.13) 

where Xi is a coordinate in plane 1, and X2 is a coordinate in plane 2, and 

{GS\T[iIjI{x2) ^2*(xi)e^"'^r(.2..i)("2+^M)'^-M-^/r(.,,,2)("2+^'') tlj,{yi)^2{y2)]\GS) (6.14) 
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were and x\ are the same coordinate as in eq (|6.13|) but hving now in plane 1 and 2 
respectively. 

Following the same steps as above, we find that the relative statistics byi is given by 



In particular, for the {m,m,n) states, the relative statistics results 6 = ^^2" „2-) ■ 

Up to this point, we have only considered the statistics for states such that the matrix k"^ 
is invertible. We consider now the case in which this matrix is not invertible. In particular, 
we study the (m, m, m) states. In order to calculate the two particle Green function for two 
particles in different layers, we have to calculate averages of the form 



I' 



Paj^e^^-e"^^-'^^^""^^-^/^'^'^'^' (6.16) 

In this expression the indices a and f3 are fixed, i.e., no sum over them is assumed. Since we 
are calculating the relative statistics, a = 1 and /5 = 2 or viceversa. The currents have the 
same meaning as in eq ( |6.11| ). In the low energy limit, we can take again the effective action 



which only contains the Chern-Simons term, and that is given by eq ( |5.1| ). In particular, 
for these states, it is more convenient to work in the basis defined by = ± a^. In this 
basis, the action is given by eq ( |5.7| ) with = 0. Therefore, the expression in eq ( |6.16[ ) 
results 

/ Va^Va- exp J ^y?[«:^-^< e,^,Q' 5+] - ^ j^^^^^ a^j^ + a;f_ (6.17) 

where k"'""'" = ~2ff(!^2m)' '^^^ integral over a+ can be performed, and gives an expression 
as the one in the r.h.s. of eq (|6.11| ) but replacing a by +. The path integral / danti+ 



represents the self-linking number of the trajectories defined by j+ . 

In the long wavelength limit, in which eq. (|6.17|) is exact, the gauge fields for the out 



of phase degrees of freedom, a~ , only enters in the linear coupling to the external currents. 
Hence, in this limit, the out of phase gauge field plays the role of a Lagrange multiplier 
field. In particular, the integral over a~ yields the constraint j~ = — = 0. Since these 

34 



currents are non-zero only on the world lines of the particles, this constraint states that the 
particles are forced to move together forming a bound state. This is a direct consequence of 
the existence of the gapless mode. If we were to keep also terms which, in the effective action 
of the gauge fields, vanish in the long wavelength limit, this conclusion will be relaxed only 
a little. By dimensional analysis, the effects of these terms should drop out like a power of a 
length scale, which typically is the cyclotron length. The only effect of these irrelevant terms 
is to give the bound states a spacial extension of the order of that length scale. In other 
words, the effect of the gapless mode in the spectrum, is to induce a long-range interaction 
that forces the quasiparticles in different layers to move together, forming a bound state. 
Notice that this is a very stringent constraint since the "constituent" "independent" fermions 
from each layer have now been eliminated from the spectrum by the existence of the gapless 
out phase gauge modes. Indeed, the fact that there is a gapless mode which propagates 
with a linear dispersion relation tells us that the first non-vanishing term in the action of 
the out of phase gauge fields has a Maxwell-like form. Thus, the Chern-Simons fermions feel 
an instantaneous force (mediated by the gapless mode) which has a logarithmic dependence 
with the distance. This is a confining force and leads to strong infrared divergencies and 
to the confinement of the Chern-Simons fermions. Only the bound states remain part of 
the physical spectrum. This picture is strongly reminiscent of the mechanism by which the 
fermionic states disappear from the spectrum of anyon superfluids p6| , ^ |. 



VII. SPIN SINGLET STATES AND SU{2) SYMMETRY 

In this section we discuss the application of Chern-Simons methods to a 2DEG which 
has an exact SU (2) symmetry. In reference a non-abelian Chern-Simons approach was 
used to construct a theory of the spin singlet (m + 1, m + 1, m) Halperin states. Here, 
we present a generalization of the approach of reference which will yield the SU{2) 
invariant hierarchies. A related, and completely equivalent approach, was developed by 
Frohlich, Kerler and Marchetti 1231 
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In the approach of Balatsky and Fradkin [|T^] (BF), instead of the two-component abelian 
Chern-Simons theory that we use in this paper, a non-abelian Chern-Simons gauge field is 
introduced. The advantage of the BF approach is that the SU (2) invariance is manifest and 
it is not the consequence of a subtle dynamical mechanism. The disadvantage of the BF 
approach is that the non-abelian Chern-Simons theory is substantially more sophisticated 
and technically more demanding than the abelian theory that we use in the rest of this 
paper. In the BF approach, the electron is viewed as a composite object which is made of 
a particle that carries the charge (the holon) and another particle that carries the spin (the 
spinon). This arbitrary separation gives rise to the existence of an abelian gauge symmetry 
(called RVB by BF). The requirement of gauge invariance forces the holons and spinous to 
be glued together in bound states, the electrons. In the BF approach, the need of a non- 
abelian gauge field is a consequence of the assignment of fractional (semion) statistics to 
both holons and spinous. In this way, SU (2) fluxes are attached to a set of charge neutral, 
spin-|, fermions, which become the spinous. The SU{2) symmetry only admits Bose, Fermi 
or semionic statistics. The holons, instead, are represented by fermions attached to U{1) 
fluxes. The U{1) and SU{2) Chern-Simons coupling constants must be chosen in such a way 
that the holons and spinous are semions. Within this approach, the FQHE of the spin-singlet 
states is the FQHE of the semions. The spin structure just sits on top of the FQHE. 

There is another, and more obvious, way to attach fluxes to particles while keeping the 
SU (2) spin symmetry untouched. Belkhir and Jain recently proposed a spin singlet 
wave function for the state with filling fraction = |, based on a composite fermion picture. 
In their construction, they attach the same number of pairs of flux quanta to both up and 
down spins (in other words, they attach fluxes to the charge and not to the spin). The same 
construction, but in the bosonic Chern-Simons language, was used before by D. H. Lee and 
C. Kane and subsequently applied by Sondhi, Karlhede, Kivelson and Rezayi in 
their theory of skyrmion states in polarized FQHE states. These approaches are represented 
in the U{1) ® U{1) theory of the preceeding sections by demanding that both up and down 
electrons see that same flux at all times. This means to choose a Chern-Simons matrix 
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which is proportional to the identity, i. e. , si = S2 and n = 0. It is easy to use this approach 
to get the SU{2) hmit of the Jain states, but with arbitrary polarization. However, it is not 
a very efficient approach to get all of the spin singlet states. The approach of BF deals with 
this states more directly. It would be desirable to construct a theory which has as particular 
cases both the SU{2) theory of BF and the U{1) theory of reference |^D[. In what follows, 
we follow the BF approach. 

More concretely, following BF, we introduce a holon field cj) and a spinon field Xa (« i) 
and represent the electron field operator ipa as ipa{x) = Xa{x). BF showed that, for an 
SU (2) invariant system, the system defined by the action of Eq. is equivalent to the 
following theory of (interacting) spinous and holons. Let the total action S be the sum of a 
charge, spin and interaction terms 

^ '^charge ~l~ '^spin ~l~ '^interaction C^'-^) 

where the action for the charge degrees of freedom is 
.Scharge = / d^x (iD'o + fx) 0(x) + j^<PKx) D^^ix)^ + J d^x ^-e ^.^a^ {x)d'' {x) 

(7.2) 

while the action for spin is 

'Sspin = jd^x [xi{x) iD^Q Xa{x) + ^xi{x) D^^Xaix)^ 

- J d'x Ae tr(6^(a;)a"6\x) + h^'{xy{x)b\x)) (7.3) 
An 3 

and an (instantaneous) pair interaction term 

'^interaction = - J d^X J d^x' ^(p(x) - p)V {x - x'){p{x') - p) (7.4) 

In equation (|7.2|) , is the statistical vector potential which turns the holons into semions. 
This condition requires that the U{1) Chern-Simons coupling constant 6 be restricted to the 
values 

1 27r 

7: = + 27r2s (7.5) 
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where semion statistics requires that m = ±2 and s is an arbitrary integer. Likewise, in 
eq. ( [7.3| ) is the SU (2) non-abehan statistical gauge field which takes values on the SU{2) 
algebra. Hence, we can expand the field in the form bf^{x) = r", where (a = 1, 2, 3) 
are the three generators of SU (2) in the spinor representation , i. e. the set of 2 x 2 Pauli 
matrices. The SU{2) Chern-Simons coupling constant k, the level of the Chern-Simons 
theory, for our system, is equal to [|l4l k = ±1. Hence, we have a ("trivial") level one Chern- 
Simons theory. All the representations of a level one Chern-Simons theory are known to be 



abelian and to correspond to abelian fractional statistics of fermions, bosons or semions |]15 
The only ambiguity left is the sign of k which is the chirality (or handedness) of the semion. 
There is a similar sign ambiguity in the coupling constant 6 in eq. ( [7.5| ). Different choices 
of these signs lead to different statistics for the bound states of holons and spinous. The 
requirement that the bound state be an electron, which is a fermion, leads to the condition 
sign(/i;) = sign(m). 

The U{1) and SU{2) charge and spin covariant derivatives are 

Dl= Id,- i{h, - I c,) (7.6) 

where A, is the external electromagnetic field and / is the 2x2 identity matrix. The gauge 
field is the U{1) gauge field which glues spins and charges together. The covariant 

derivatives have been chosen in such a way that holons and spinous have opposite charge 
with respect to c^. Hence, the strong fiuctuations of this field binds holons and spinous into 
states which are locally singlets under the "RVB" gauge transformations, i. e. , electrons. In 
fact, all the gauge field does is to enforce the constraint that the 3-current of the holons 
equals the 3-current of the spinous, as an operator statement in the physical Hilbert space. 
This is seen clearly from the equation of motion generated by 

J^^^(x) = /f- = =^ 0t(^)0(^) _ ^ xl{x)Xaix) = (7.7) 
Sco{x) ^t^^ 

Thus Ne, the number of charges, and and Ni, the number of up and down spins, must 
obey the obvious relation A^^e = + N^. Also, eq. (|7.7|) tells us that the local particle 
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density operator p{x) can be identified with the holon charge density operator, i. e. , p(x) = 

We will not attempt to go into the details of the full non-abelian theory. Rather, we 
will use it to determine the allowed fractions for SU{2) invariant states. Thus, we will just 
consider the AFA equations for this theory. There are two sets of AFA equations, one for the 
charge sector and one for the spin sector. The AFA equations for the charge sector are just 
the AFA equations for a charged interacting semion liquid with A^^e particles in an uniform 
magnetic field with flux quanta, z. e. , a FQHE of semions. A simple application of the 
methods of reference yields the constraint for the charge density operator 

6S 



6ari(x) 



=^ 3o{x) = -eB,{x) (7.8) 



where joi^x) is the charge density operator. This equation, when specialized on fluid states, 
for which p and Bc{x) are the average density and average U{1) statistical charge flux 
respectively, yields the allowed filling fractions. Similar considerations for the SU{2) gauge 
field 6° (with a = 1, 2, 3) yield a constraint for the local spin density operator 

=^ j-(x) = V^t^V. = (7.9) 



where is the SU (2) spin flux. 

For fluid states, we generalize the Average Field Approximation and replace these exact 
local operator identities by translationally invariant averages. This replacement may be 
problematic in the case of the RVB gauge fleld since its only mission is to enforce exactly 
the constraint eq. ( |7.7| ). At the level of the wave functions for the allowed ground states, 
this constraint simply means that every coordinate for a charge degree of freedom has to 
coincide with the coordinate of a spin degree of freedom. 

Thus, we seek fluid states with N^. particles, flux quanta and flUing fraction v = 
Ne/N^. We wish to determine the filling fractions and spin for which the ground state is 
a fluid. In addition to the usual AFA equation for the charge, that will give the allowed 
fractions, we will now get conditions for the spin and polarization of the allowed states. 
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In the charge sector, we have a FQHE of spinless semions. In this case, the AFA consists 
of a system of A^e spinless fermions filUng up effective Landau levels, exactly as in our earlier 
work for spin polarized electrons (and anyons!) 0. The same line of argument that was used 
to derive the allowed fractions in Section |11 A| now tells us that, from eq. ( |7.8| ), the effective 



number of fluxes is 



= (7.10) 



and, hence, the allowed fractions, u^, for the SU{2) fluid states satisfy 



where p is a positive integer and the + sign corresponds to a "particle-like" FQHE (^. e. , 
the effective flux parallel to the external flux) while the — sign holds for a "hole-like" FQHE 
{i. e. , the effective flux anti-parallel to the external flux) . By using the allowed values of 
9 we find that the allowed fractions for the SU (2) fluid states are of the form 



±f \ mp _ 2p 
u [p, s;m) = = — -— ('-12) 

[2sm — l)p ± m ±2 + (4s — sign(m))j9 

where we have specialized for the case of interest, m = 2sign(m). 

The hierarchy of FQHE states of eq. ( [7. 121) is a generalization of the states found by BF, 



which are obtained by setting p = +1 and m = +2, i. e. , z/+(l, s; —2) = = 2, |, |, . . . , 
which coincide with the Halperin-Haldane spin singlet states. 

The state with z/ = | is also part of the hierarchy of eq. ( [7.12|) , where it is realized 



as the state with p = 2 (—2) for m = +2(— 2) and s = 1. Numerical studies show that, 
for Coulomb-like interactions, this state is not favored and that the (3, 3, 1) is an accurate 
representation of the ground state. 

Finally, the this hierarchy has a state at filling fraction |. This state is found as the p = 
— 10 ("hole-like"), m = —2, s = 0, member of the SU{2) hierarchy, or as the pi = P2 = —5, 
Si = S2 = s and n + 2s = 1 member of the U{1) U{1) hierarchy. This is the only spin 
singlet state that has yet been seen experimentally. 
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In contrast with the states that we found in section ^IA\ using the abehan theory for 
bilayers, in the SU (2) theory, the states at each filling fraction, are arranged in irreducible 
representations (multiplets) of SU{2). Thus, the SU{2) theory does not have any redundant 
states. The spin and polarization of the states is determined by eq. (|7.9|). The only subtlety 
here is that, since the components of the total spin do not commute with each other, one 
can only determine the total spin and total projection. For a system without a boundary, 
the choice of total spin S and total projection along an arbitrary polarization axis, say 
Sz, are constants of motion which are invariant under local SU{2) gauge transformations 
(but, of course, change under global SU{2) rotations). There are two generic situations of 
physical interest: (a) spin singlet states (or with microscopic total spin S/N^ ^ 0{1/Ne)) 
and (b) states with macroscopic spin, S ~ A^e, i- e. ferromagnetic states. Thus, the total 
z-component of the spin polarization M = |(A^| — Ni) obeys 

2M = N^-N^ = ^{Bt)L' (7.13) 

where is the area. It is clear that it is possible to construct all multiplets with spin 
\S\ < In the thermodynamic limit, the spin singlet states have S = and, hence, 
N-^ = Ni = Ne/2. In contrast, the ferromagnetic states have, with an appropriate choice of 
the quantization axis, a non-vanishing extensive value of M and, hence, a non-zero value of 

The wave functions for the spin sector of the spin singlet states have to be determined 
from the states of a level one SU{2) Chern-Simons gauge theory with sources in the 
fundamental representation. It was shown by Witten [|15| that these wave functions are 



correlation functions of conformal blocks of a conformal field theory in two Euclidean di- 
mensions, the SU (2) level one Wess-Zumino- Witten model. This fact was used by Read and 
Moore |T6| and by Balatsky and Fradkin |1^ to show that, the wave function of the spin 
singlet FQHE has a factor which is precisely this conformal block correlation function. It 
was also noticed that this factor coincides with the Kalmeyer-Laughlin wave function 



3^ for a Spin Liquid. 
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The states with macroscopic spin have a somewhat different physics. The existence of 
a non-zero average field should make a mean field approach more sound. The spin sector 
of this mean field theory has N-^^i = ^Ne ± M spin up and spin down spinous each feeling 
an effective uniform magnetic field of ±^M. Notice that, because of the SU{2) invariance, 
there is no Zeeman term and only the orbital degrees of freedom see this spin-dependent 
external field. It is easy to see that the highest weight ferromagnetic state with maximal 
spin is obtained by filling up the lowest Landau level of the up spins while leaving the down 
spin sector empty. 

The charge and spin sectors are not decoupled from each other. Firstly, the constraint 
of eq. ( [7.7|) sets the local charge density to be the same as the local spin density. The wave 
functions of the allowed states have to satisfy this local property. Secondly, if the system 
is SU{2) invariant, all of the states in a given SU{2) multiplet must have the same filling 
fraction. Since the fully polarized states have to span all of the Jain states for a single 
layer system, we must conclude that SU{2) states which are not in a main Jain hierarchy 
cannot achieve the maximum polarization. In other terms, there is an upper bound for the 
spin polarization and, hence, for the total spin itself. Thus, the filling fraction u and the 
spin S of the state cannot be set completely independent from each other for the allowed 
states. In other terms, there should exist a set of selection rules which determine the allowed 
combinations of total spin and filling fraction. Similarly, it should be possible to construct 
a unified theory of all the FHQE states with SU (2) symmetries, instead of the apparently 
separate descriptions of the spin singlet and the fully polarizable states that we use here. 
We will return to these issues in a separate publication. 

We conclude this section with a comparison of the states that are obtained by this SU{2)- 
symmetric approach and the U{1)®U (1) theory that we use in the rest of this paper. A direct 
inspection of the allowed fractions eq. ( p.l3| ) and eq. ( |7.12| ) for the U{1) ® U{1) and SU{2) 
theories respectively, shows that they do not yield the same allowed fractions. For instance, 
the "Fermi Liquid" (compressible) states, with the same occupancy of the two layers, allowed 
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by eq. (|2.13|) have filling fractions ^ = 2, 1, |, ^, |, |, |, i, |, . . . (r = 1, 2, . . .). In contrast, the 
allowed SU{2) "Fermi Liquid" (compressible) states are = 2, |, |, |, f , • • •• Clearly, the 
fractions ^ (with A; = 1, 2, 3, . . .) cannot be realized as SU{2) compressible states. Among 
the states which appear in both hierarchies, we find an incompressible spin unpolarized 
state at filling fraction ^. It is worth noting, that there is experimental evidence of an 
incompressible state at ^. It is not possible to construct a fully polarized Jain state with 
this filling fraction although it may be constructed as a hierarchical state. It is strange 
that this is the only observed fraction for which a hierarchical construction is needed. It is 
believed that the experimentally observed state is polarized. 

We have further checked that all of the states in levels 1 and 2 in the SU{2) hierarchy 
span the entire level 1 t/(l) ® U{1) states. Similarly, the level 4 SU{2) states span the level 
2 U{1)®U{1) states, and the levels 3 and 6 SU{2) states span the level 3 U{1)®U{1) states. 
However, a large number of incompressible f/(l) ® U{1) states cannot be realized as SU{2) 
states. In a way, this should not be surprising since the U{1) ® U{1) symmetric theory 
may only generate SU (2) as a dynamical symmetry. Nevertheless, it is a puzzling result 
since, in this discussion, the form of the interaction terms has not entered and, hence, the 
symmetries of the Hamiltonian have not had a chance to play any role yet. However, even 
more surprising is the fact that not all of the SU{2) states can be realized as U{1) ® U{1) 
states. It is easy to check, for instance, that the allowed SU{2) state with = y has no 
counterpart in the U{1) ® U{1) states (unless polarized states are also considered). 

VIII. CONCLUSIONS 

In this work, we generalized the fermionic Chern-Simons theory for the Fractional Hall 
Effect (FQHE) which we developed before, to the case of double layer-systems. 

We studied the semiclassical approximation around the liquid-like mean field solutions. 
We found that we can describe a hierarchy of double-layer FQHE states which include 
the (m, m, n) states, and the (m, m, m) states, with filhng fractions u = and z/ = — 

\ 1 1 / 1 \ 1 1 J 1 o m+n m 
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respectively. For all these states, the mean field ground state has a gap to all the excitations 
and the semiclassical approximation is a well controlled perturbative expansion. Within the 
solutions of the saddle point equations, we also encounter the generalization to the double- 
layer systems, of the compressible states discussed by Halperin et al [|1^]. It turns out that 
the liquid-like solutions of the saddle point equations predict a very rich phase diagram for 
double-layer systems. The structure of this phase diagram should depend on the particular 
microscopic form of the pair interactions. 

We studied the electromagnetic response functions of these systems and calculated the 
experimentally accessible optical properties. In general, for the so called (m, m, n) states, 
we found that the spectrum of collective excitations has a gap for all the modes. We 
derived the absolute value squared of the ground state wave function. We found that it has 
the Jastrow- Slater form, with the exponents determined by the coefficients m, and n. We 
also found that the (m, m, m) states, have a gapless mode which may be related with the 
appearance of interlayer coherence [0,|lT|. Our results also indicate that the gapless mode 
makes a contribution to the wave function of the (m, m, m) states analogous to the phonon 
contribution to the wave function of superfluid He4. This factor is crucial to obtain the 
correct spatial correlations from the ground state wave function. 

In all the cases we verified that the density correlation function saturates the /-sum rules 
associated with the conservation of the number of particles in each layer separately. 

We calculated the Hall conductance and verified that it gives the correct value already 
at the semiclassical level of our approach. 

We calculated the charge and statistics of the quasiparticles. We saw that the charge 
is determined by the filling fraction in the layer divided by the corresponding number of 
effective Landau levels filled. We found that the statistics of the quasiparticles is well defined 
only for states which satisfy that the Chern-Simons matrix of the effective action for the 
statistical gauge fields is invertible. For instance, for the (m, m, n) states, the statistics 
for a quasiparticle in a given plane is determined by ■T^ ^ , and the relative statistics is given 
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by o" 2 - On the other hand, for states such that is not invertible, as for instance the 
{m, m, m) states, we found that the statistics of the quasiparticles is not well defined. In 
particular, we saw that the effect of the gapless mode is to induce a long-range attraction 
that forces the quasiparticles in different layers to move together, forming a bound state. 

We have also compared the results of this U{1) (8) theory with a generalization 
of the SU (2) theory of Balatsky and Fradkin. We showed that, for many fractions, the 
C/(l) U{1) can be embedded inside an SU{2) theory. In this way we were able to identify 
SU (2) states in our theory for bilayers. However, we also found that the two hierarchies are 
not completely equivalent. 

Finally, we would like to remark that, within our approach, the gap of the collective 
modes is always proportional to the effective cyclotron frequency. In particular, this holds 
also for the out of phase modes. In principle, one expects that the zero momentum frequency 
of these modes should be determined essentially by the interlayer pair potential, V12, but, as 
in the single layer problem, we can not obtain this result within the gaussian approximation 
and we expect the non-gaussian corrections to give a substantial correction to this energy 
gap as well. 
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APPENDIX A: 



In this appendix we follow the method described in Appendix A of reference to prove 
that both theories yield the same physical amplitudes (i.e., that they are equivalent). The 
idea is to compute the same (arbitrary) amplitude in both schemes. We do the calculation 
in the path-integral language. We can follow step by step what we did in reference [^. The 
only difference is that in this case, there are two species of fermions. The fermions living 
in one layer are identical, but they are distinguishable from the fermions living in the other 
layer. 

The standard action for nonrelativistic particles coupled to the electromagnetic field is 

/ + 00 f 
dt dhj^{z,t)A>^^ + S^^tter (Al) 

where /i = 0, 1, 2 is the space-time index, a = 1, 2 labels the layer, and ^matter includes both 
the kinetic energy of the particles as well as their pair interactions. The boundary conditions 
are 

lim Zjit) = Xi 

where the final states are a permutation of the initial states. The current j^(^, t) is a three- 
vector of unit length tangent to the world lines and takes a non-zero value only on the world 
lines of the particles (r(a)) 

dz- 

r{z,t) = Y.mt)-m)^ m 

If we couple the above system to the Chern-Simons fields a^, the action becomes 

5new = S + Y^ r^dt f dh j^{z, t)a'^ + Y.^ f d^z e^,A<9a^ (A4) 
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It is possible to perform the functional integral over the statistical gauge fields exactly. These 
fields enter in only two terms of S^ew'- the current term and the Chern-Simons term. Thus, 
the average over all the configurations of the statistical gauge fields has the form 

(exp(^ Y: I dhr^{z)a^^{z)))cs ^ exp{tl[f^]) (A5) 

where the notation (0)cs indicates the average of the operator O over the statistical gauge 
fields with only a Chern-Simons action. 

Our goal is to find out which has to be the form of the Chern-Simons coefficient, i.e., of 
the matrix Kq^, for the r.h.s. of eq (|A5|) to be equal to one. 

After integrating out the Chern-Simons fields we obtain the following expression for I[j^] 

nf,] = E / d'z J d'z'r^{z)G>^'^{z, z')fM') (A6) 

a/3 

where G^y{z, z') is the Green function 

G^,{z,z')=e,,xdt^Go{z,z') (A7) 
and Go is the Coulomb green function 

-d^GQ{z,z') = 6^''\z- z') (A8) 



By direct substitution of eq ( |A8| ) into eq ( [A7| ), we can write the exponent /[j^] in the form 



m = E jd^j dz r^iz)e^'''di'^Goiz, z') /^{z') (A9) 

0/3 

Using the magnetostatic analogy, we can now regard j° as a current in three dimensional 
Euclidean space and use it to evaluate the expressions in eq ( [A9| ) . Let C° be a vector field 
related to by the equation ("Ampere's Law") 

VxC'^=r (AlO) 

such that 

V ■ = (All) 
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The solution of eq ( |A10| ), subject to the constraint eq ([AllD , can be found in the same way 
as in reference M. The field is given by 



(A12) 



By substituting eq ( |A12D back into eq ([A 91), we find that /[j^] takes the simpler form 



E 

a/3 



2 , d^zr^{z)3;{z) 



(A13) 



Now, since the currents are non-zero only on the world lines, we can rewrite the volume 
integral in eq ([A13|) in the form of a line integral over the configuration F". The set of closed 
loops are the boundary a surface S", F" = We can then apply Stokes' theorem to 
get the result 

-l\Q/3 



m] = E 



a/3 



n'^ jjl 



(A14) 



where is a vector field normal to the surface S". The integral /^c da n° is an integer 
which counts the number of times the current pierces the surface S". 

The result of eq (|A14| ) means that the average over the statistical gauge fields eq (|A5|) 
has the simple form 

(exp(z Y: I d'zr,iz)A^,iz)))cs = exp ^ ^^^^ / da f^) (A15) 

We are now in position to ask when are the two systems equivalent. In other words, when is 
the average over the statistical fields equal to one or, at least, independent of the integral in 
eq ( [A14| ). In the first case we would have proven that the amplitude in the system with the 
statistical gauge fields is exactly equal to the same amplitude calculated in their absence, 
whereas in the second case all amplitudes differ by a constant phase factor, which can be 
dropped. By inspecting eq (|A15D we see that the amplitude is equal to one if 



27r 



^ 2si 



V 



n 2s2 



(A16) 



where Si,S2, and n are arbitrary integers. 
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APPENDIX B: 



In this Appendix, we derive an expression for the energy of the ground state which 
includes the corrections due to the gaussian fluctuations of the Chern-Simons gauge fields. 
We evaluate this energy for two different incompressible states with filling fraction v = l.We 
show that the mean field degeneracy is lifted. In particular, within our approximation, we 
show that the (1, 1, 1) state is the lowest energy state. 

We begin with the derivation of an expression for the ground state energy. The quantum 
partition function for this problem at zero temperature is given by eq ([2.81). In order to 



calculate the ground state energy, we will set the external electromagnetic field to zero. Since 
the action is quadratic in the fermions, they can be integrated out. After expanding the 
Chern-Simons gauge fields around their mean field solution < a° >, the partition function 
can be written as follows 

Z = e'^^^ J Va^ e^^^ff(^2) (Bl) 

where S'mf coincides with the expression for the effective action in eq ( p.9D but setting 
a'^ =< > and = 0, and Scs{(i'^) is given by eq (|2.17|) with = 0. If we define 
II^^(x,?/) as the sum of the polarization tensor n^^(x,y) plus the tensors corresponding to 
the interaction term and the Chern-Simons term in eq ( |2.17| ), we can write the effective 
action in eq ( plD as 



^eff(a^) = lld'xl d'y a;(x) ^^^{x, y) ~al{y) (B2) 

Since this action is quadratic in the fluctuations of the Chern-Simons gauge fields, they can 
be integrated out explicitly. The partition function becomes 

where Vetijl^p) is the functional determinant of the gaussian fluctuation operator. Thus, 
we can write Z = e*'^*°* where 



49 



•S'tot = <S'mf + - lnT'et(n: 



2 a/3J 



Smf + -Trln(n 



af3) 



i f (PQ f duj 



Smf + 7; TTTT^ / TTTT 1^ det(n^^) (B4) 



In eq. ( P4|) TiA stands for the functional trace of an operator A while trM is the algebraic 
trace of a finite matrix M over the indices a, (3, fi and u. Similarly, det M is the algebraic 
determinant of the matrix M. 

In the limit of infinite time span, T ^ 00 (or zero temperature), the quantum partition 
function can be written in terms of the ground state energy Eqs, a.s Z = e~'^^^^'^^, where 

is the area of the system. Therefore, using eq (p4|) , we find that the ground state energy 
is given by 

Scs^^^'f- 5^/11/ (^ln(datn-) (B5) 

It can be shown that det II[^^ is equal to the square of the denominator of the density 

^ - r ^ -|2 

correlation function, D{uj,Q), i.e., det 11^^ = D{uj,Q) . Here, as in the single layer case, 
the zeroes of the denominator of the density correlation function determine their poles. 
Therefore, the zeroes of D{uj, Q) are the collective modes of the system. In other words, we 
can write D{u!, Q) = Yixi^'^ ~ ^x), where are the collective modes of the system. 

To evaluate the integral in the r.h.s of eq ( p5|) we consider a cutoff function f/<^{u!,Q) 
such that fhioj^Q) — > when u — > ±00, and d^fiyioj^Q) 0,2. e., it is an adiabatic cutoff. 



Therefore, the integral in the r.h.s of eq (B5) can be calculated as follows 
/ |^lndet(nr,) = 2/ ^f^{u,Q)\nD{u,Q) 



1 7^^^(^'^)7;7^^-^(^'^) 



777;t/a(w,<5)— — 

(27r) D{uj,Q) 

where the last equality was obtained integrating by parts. Using this result eq ( [B5D becomes 

Eos - E^i ^ ^ / H / '3)^«-«(-. Q) m 
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After integrating in uo we get 

^GS - ^GS^ = / 1^ E -a/a(Q) (B8) 
In eq. (|B8| ) tlie sum over modes includes the case of eventual multiplicities and degeneracies. 



We will evaluate this expression for the state I, which has u = 1 and pi = P2 = 1, 
■Si = S2 = and n = 1, and for the state II, which has also u = 1 and pi = P2 = 1, but 
Si = S2 = 1 and n = —1. It can be shown that the collective excitations of the state II 
coincide with the ones of the state (3,3, 1) (derived in Section |III A| ), but replacing vu by 



— fi2, and using that uj = ^ and B = ^. Substituting the expressions for the collective 
modes calculated in Section |III B| for the state I, and the ones calculated in Section [III A| for 



the state II, (but replacing V12 by — fi2, w = ^ and B = y), into eq (p8|), we obtain 

ELs - Eg = -^Yl^ I dx[-A - (1 + M)- + {2x{l + M)r^]U{x) (B9) 

where M = ^{vn — ^12), x = and we have kept only linear terms in x. In particular, 
using a sharp cutoff to calculate the integral in the r.h.s of eq (|B^), such as < 2B, we 
find that Eq^ — Eq^ < 0, in both cases, for Coulomb and for short range interactions. 

A few comments on this result are in order. Firstly, since this theory does not have a 
small parameter, it may be argued that a result obtained in the leading, gaussian, correction 
to the mean field result might be affected significantly by higher order corrections. For a 
generic pair interaction, this should be the case and the leading order result may not be so 
significant. Notice, however, that the approximations that we have made are accurate only 
for interactions whose characteristic length scales are smaller than the cyclotron radius. In 
this situation, the gaussian effects should become dominant and the selected state should be 
unique. However, if the interaction is important only at distances long compared with the 
cyclotron radius, the selection of the states may be different and a complex phase diagram 
should then be expected. But, in such cases, other states , such as crystals, may also compete 
effectively with the fluid states. 

Thus, we have shown that, for an interaction with a sufficiently short range, the state 
(1, 1, 1) is the fluid ground state. We believe that similar arguments can be made in all the 
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other cases where we found multiple solutions (with the same caveats about the range of 
the interactions). 
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